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1 Introduction 



Matrix models have been studied intensely since the classic paper pQ, where it was 
realized that they enumerate planar diagrams. Remarkably, it was recently discovered 
that there is a description of Chern-Simons (CS) theory 013] and Holomorphic Chern- 
Simons (HCS) theory on certain manifolds as particular matrix models. Whilst 

interesting in its own right, this discovery also has deep implications for supersymmetric 
gauge theory in four dimensions. 

The connections between CS/HCS theory and field theory in four dimensions is the 
following (see [7] for a review). Type II string theory on a Calabi-Yau (CY) manifold 
can be twisted in one of two ways to give a topological string, the so called A and B 
models [Hj. Furthermore, it was shown in [H] that the open string A-model on the CY 
manifold T*M, where M is a Lagrangian submanifold, is equivalent to CS theory on M 
with gauge group determined by the Chan-Paton factors. It was also shown that the 
open string B-model on a CY is equivalent to HCS theory on that CY. So calculating 
the partition function of these CS/HCS theories amounts to solving the topological 
A/B string. 

For certain topological correlators the A/B twist is trivial and thus these topological 
correlators give physical correlators, this in fact happens for precisely the correlators 
that correspond to F-terms in the resulting four dimensional gauge theory |Sj. So the 
conclusion is that the superpotential can be computed entirely from the topological 
string. This was made extremely precise in [TUIId] , where it was shown that the effective 
superpotential in four dimensions W e ff is given by, 

W eff = N^-rS, (1.1) 

where, N is the amount of RR flux, F is the leading order contribution to the topo- 
logical string free energy, S is the gaugino condensate and r is the gauge coupling. So 
showing that topological strings can be described by matrix models is very suggestive 
that the matrix model structure can be uncovered directly in field theory. For the case 
of SU(N), M = 1 gauge theory with adjoint matter, this was found in [TH IT2]. 

After the initial work of Dijkgraaf and Vafa, it was shown how to generalize their 
work to the other classical gauge groups from several points of view. Matrix model 
generalizations were considered in the perturbative supergraph techniques of 

[TT] were considered in [T3] and generalized Konishi anomaly techniques of [T2j were 
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considered in [TniUZj- All these works were studying J\f = 1 SO/Sp gauge theories in 
four dimensions with adjoint matter and single trace superpotential. 

The present work is concerned with how the calculations of jTHI (which will be 
reviewed in section 2) can be performed in the CS model matrix model of j2HHJ- The 
main result of this paper is the calculation of Fi in the CS theory on S 3 with gauge 
group SO/Sp. The partition function of this CS theory has in fact been calculated to 
all orders in [IE], the present work explores by explicit calculation, the matrix model 
description of CS theory. At first glance there appears to be an inherent contradiction 
between a naive extrapolation of the results of and the known partition function 
|18j . We will find in this paper that although the CS case is more complicated than 
the Lie algebra matrix model, the matrix model realizations of CS is not incorrect. 

In P], it was shown that the CS model matrix model is of the same type as the 
B-model matrix model but with a rather complicated double trace potential. It must 
be because of the double trace that the B-model calculation does not translate to the 
CS model. 

This paper is organized as follows. In section 2 we review the calculation of ^3] 
and present a new way to obtain the same results. This new method will generalize 
to the CS matrix model. In section 3 we describe the topological string and Calabi- 
Yau geometry which is inherently being studied. Section 4 will contain a discussion 
of the free energy of SO / Sp CS theory, what the naive contradiction is and what we 
will calculate from the matrix model. Section 5 will contain a derivation of the loop 
equation for matrix models with Haar measure, something which has not appeared in 
the literature. We will see explicitly why the method of [T3j breaks down. In section 
6 we will calculate the leading order free energy of the CS matrix model for groups 
SO/Sp, this is a straightforward generalization of j3] but is included for completeness. 
In section 7 we will calculate the subleading free energy, this is the main result of 
the paper. In section 8 we will discuss the four dimensional gauge theory which this 
analysis actually corresponds to. 

2 Matrix Models for Classical Groups 

Matrix models for all the classical groups were first considered in |19j, where they 
wrote down the appropriate measures in eigenvalue form. More recently, the authors 
of ^3] studied the orientifolded CY geometry of [I] and the resulting matrix model. 
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As is well known, orientifolding internal geometries in string theory leads to the gauge 
groups SO/Sp. It was shown in |i3. that the relevant matrix model, when in eigenvalue 
form, has a measure on the Lie algebra of 5*0 / Sp, as considered previously in ^H] • 

2.1 A first look at F x 

Techniques for calculating 0(N~ 2 ) and higher corrections to the Hermitian matrix 
model free energy were considered previously in j27Hl21|(thev correspond to four di- 
mensional gravitational F-terms). The same method was used in to derive the 
0(N~ l ) correction to the free energy (Fi) and will now be reviewed. First, lets set up 
some notation. The partition function of the matrix model is, 




(2.1) 



where, W = YujLi §^ 2j > an d $ is in the adjoint of SO/Sp. For a single cut model, 
the number of eigenvalues in that cut is M, and we define S = g s ^- The resolvent is 
defined as, (g is for genus, c is for crosscap) 




(2.2) 



Now we need the relation, 



where -4t(x) 




(2.3) 




(2.4) 



to see that 




(2.5) 



Then we derive a loop equation for the full resolvent, the answer being 



2 



£ 



dx W'(x) 





^(^) + g 2 s — (x)uj(x), 



(2.6) 



27ri x — x 



x 
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and using the g s expansion of uj(x), we extract the zeroth and first order loop equations, 



/ dx' W'(x') . . , . , 1 . . 

2 & — ^-Wx' = 2uj 1 (x)lu (x) - -ooo(x). 2.8 

J c 2m x — x x 

Now one observes that a solution for uj\ (x) is given by, 

U1{X) = Yx~A-dS (2 - 9) 
then using (|2.5jl we see that this implies 

F, - -if. (2,0) 

And so we have the first correction to the free energy in terms of the planar free energy, 
valid for any single trace potential. 

2.2 A second look at F 1 

The above method of obtaining F\ does not generalize to the case of CS matrix models, 
as will become evident in section 5. So we will now find F\ for this model using a rather 
different method than the one above, one which will generalize to the case of the CS 
matrix model. 

We will assume there is just a single cut (—a, a). We then modify the potential as 
follows, 

W(x) = W(x) - yln(x). (2.11) 

We will denote the free energy of this modified model F. This modification effects the 
loop equation (|2.6|) . by adding to the LHS, a term 

9s i dx' uj{x') ^ 



2 J c 2ni x'{x — x') ' 

where as usual, the contour does not encircle the point x. Now by deforming the 
contour to infinity we pick up the point x but the integral around infinity vanishes and 
we see that 

U{X ' ] =gM. (2.13) 
J c 2m' x'{x — x') x 
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So modifying the potential in this way leaves the leading order loop equation ()2.7|) 
unchanged, but the 0(g s ) loop equation (j2.8J) becomes, 

dx' W'(x') , , . , . j k 

— — ^-f wi x = W! x Vo * 2. 14 

c 2m x — x 

which has only the solution io\ = 0. This implies that F\ is a constant which can be 
taken to be 0. 

The next step is to relate the free energy of the modified matrix model F to the 
free energy of the unmodified, CS model F. Since the leading order loop equation is 
unchanged, F = F . Now we will need to introduce the density function p(x), it is 
related to the resolvent by 

p(X) 



uj(x) = / ^-^-d\. (2.15) 
x — A 



By the saddle point approximation, we can see that 



Fi = ~ / p(X)\n\d\ - F l (2.16) 

* J— a 

i r r A i 



, ; P(A) (P / -—^dx)dX-F 1 , (2.17) 



1 

2 

1«9F 



o 

A 

wb(x)dx-Fi, (2.18) 



-F v (2.19) 
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where A is some large cutoff. Casting F\ as an integral of a 1-form over a Riemann 
surface is natural in the context of matrix models with flavour which also has an 0(g s ) 
correction to the free energy j22]- So (|2.19J) implies that 

« - "if- < 2 ' 20 > 

agreeing with the previous derivation of the same result. Generalizing this procedure 
to the CS matrix model will be the focus of this paper. 

3 Topological String Geometry 

Chern-Simons theory on S 3 with the gauge group SU(M) describes topological A- 
branes wrapped around the S 3 in the deformed conifold T*S 3 . After the appropriate 
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involution of the S 3 conifold geometry [TH] the gauge group of the CS theory gets 
replaced by SO(2M) or Sp(2M) depending on the sign of the crosscap. The involution 
goes through the usual web of dualities and large N transitions. After the large N 
transition the closed string geometry is an orientifold of O(-l) + O(-l) -> P 1 , the 
IP 1 becoming an HIP 2 . It is important that the involution does not have any fixed 
points. The mirror B-model geometry is again an orientifold of some deformed CY. 
The involution now has fixed points, two orientifold planes. This geometry has been 
also considered in [23] ■ It can be viewed as the end point of the large N transition 
on the B-model side when the B-model branes that are mirror to the A-model branes 
on S 3 disappear leaving only two orientifold planes. This two orientifold planes give 
a subleading contribution F\ to the free energy that is not present for SU(M) gauge 
group. F\ counts the holomorphic maps of HP 2 into the resolved conifold. 

After the canonical quantization the CS theory can be reduced to certain matrix 
model integrals [2]. Unlike usual Hermitian matrix models where the integration is 
performed over the Lie algebra measure, the integration in the CS matrix model is 
over the Lie group. The matrix model also can be viewed as a result of canonical 
quantization of the HCS on the B-model side with a potential that contains double 
trace terms [3]. The mirror of T*S 3 is given by the blownup of 

xz= {e u -l){e v -1). (3.1) 

Clearly, the imaginary u direction is compact, with period 2iri. The appearance of the 
group measure in the matrix integral can be interpreted as a result of the counting of 
all images of the D-branes in a matrix model with the Lie algebra measure [3] 

nn^ - u i + 2?rm ) 2 ~ ri sinh2 ( !fi T^)- ( 3 - 2 ) 

n i<j i<j 

If the gauge group is SO(2M), each eigenvalue Ui of the matrix has its partner —u^, so 
an additional product with the plus sign between eigenvalues appears in the measure. 

It is possible to take the planar limit for this matrix model and obtain a spectral 
curve, in fact the spectral curve is the nontrivial part of the B-model geometry which 
is obtained after the large N transition. The leading contribution F to the free energy 
is obtained from the integral of the resolvent over a noncompact (B) cycle. It is almost 
trivial to generalize this to SO(2M) matrix model. The main concern of this paper 
will be the subleading part of the free energy (i 7 i),that is due the orientifold planes 
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and is absent for SU(M) case. From the closed string point of view, the answer must 
be an integral of a meromorphic one-form from the orientifold planes to some fixed 
point at infinity. From matrix model point of view this one-form must be related to 
the resolvent. 



4 Free Energy. 

In this section we review known results for the free energy and relate the parameters of 
topological string and the matrix model. The A-model orientifold of the conifold T*S 3 
has been considered before. In [TSj the partition function of SO/Sp Chern-Simons on 
S 3 was calculated to all orders in g s and given a closed string interpretation. There it 
was found that F^° /Sp = 1/2F 5C/ , i.e. 

fi so/ *« = i£^. (4.D 

n=l 

The coefficient of \ is consistent with the orientifold action. From a closed string 
perspective, F 1 includes only the holomorphic maps of HP 2 into the resolved conifold 
which are odd wrappings, i.e. only Z2 equivariant maps contribute to the instanton 
expansion. Furthermore the area of an HIP 2 instanton is half that of a S 2 instanton. 
So 

F so,s P{t) = ± £ (4 2) 

n odd 

where the +(— ) sign is for SO(Sp) respectively. In the SO(2M) and Sp(2M) matrix 
model, we will use a t 'Hooft parameter S = g s M, related to the Kahler modulus t, 
the size of blown-up HP 2 , by 

t = 2S±g s . (4.3) 

This implies that the following relationship between Chern-Simon's free energy F cs 
and the matrix model free energy F MM , 

F MM (S) = F£ S (2S) (4.4) 
F MM {s) = F cs {2S) ^_dF^Sl (45) 

So more explicitly we have, 

F 1 MM (s) = ±\-y:— + > ~i ac» 
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with +(— ) sign for SO(Sp) respectively. So we immediately see that the derivative 
relation (j2.10j) does not hold for the A-model, this was the observation which motivated 
the present work. Note that the matrix model subleading free energy has two pieces. 
One comes from the nontrivial relation between t 'Hooft parameter S and the Kahler 
modulus t, the other is the contribution of the orientifold planes. 

Now another scenario where 0{g\) corrections appear is gauge theories with funda- 
mental matter. In j22j, the authors found that F\ is an integral over the spectral curve 
of the leading order resolvent ujq, consistent with what one would expect from [SUES]- 
The fundamental matter shows itself in the matrix model as a subleading term in the 
potential. We will see that the contribution from the orientifold planes also comes in 
as a subleading term in the matrix model potential and therefore also can be expressed 
as an integral of the resolvent. 

5 Loop Equation 

As in [T3j, we derive a loop equation needed to find the leading and subleading order 
resolvents, ujq{z) and u)\(z). Since the matrix integral is over the Lie algebra group 
rather than over the Lie algebra the measure factor is different, and as a consequence 
the expression for the resolvent in terms of eigenvalues is different. This does not 
change much in the derivative relation between ujq(z) and U\{z). This relationship 
does not however, appear to lift to a nice relationship between F Q and F% as ([2.100 . 



As in the matrix model with measure on the Lie algebra, the integral of the resolvent 
must be compatible with the log of the measure so we define, 



Z 




-)sinh 

2 ; 



2 ,Ui + Uj 

[ 2 



)exp 





(5.2) 



When the group is SO(2M), this becomes, 



UJ 




(5.3) 



It behaves as 



uj(x) = ±2S + 0(0 



x 



±oo. 



(5.4) 
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We restrict ourselves to one cut solutions so, we assume that u(z) has one cut that 
runs from —a to a along the real axis. We have defined S = g s M. With a potential 
given by W(x) = Y^jLi ^f x ^ ^he relationship between the resolvent and the free energy 
F is 

^(x)=2Scoth(|) + ^-yF, (5.5) 

where the differential operator can be worked out by Taylor expanding coth (f-^-) 
around x. The resolvent and the free energy have expansions in g s given by, 

g,c 

co(x) = J29 2 s 9+c ^ 9+c (x). (5.6) 
Combining (|5.5jl and ()5.6|) . we see that, 



-o = ^o + 2Scoth(f), 



d 
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' ~ dV(x 



■F,, j>0. (5.7) 



The loop equation for this model can be derived by demanding reparametrisation 
invariance of the partition function (details are in the appendix). It is given by, 



\j\x) - g s coth(:rV(x) + 2g s S - 2S 2 - 1Cu;{x) + 9 l^—u{x) = 0. (5.8) 
Where JC acts as, 

Kf(x) = j c ^~ coth(^)W'(z)/(z). (5.9) 

The contour C encircles the cut but not the point x. When we insert the expansion 
()5.6|) into the loop equation (|5.8|) . the first two equations we get are, 



O(g s ) : -c 2 (x)-2S 2 = ^ (x), 

0(g\) : Ui(x)uj (x) + 2S — coth(x)u;o(x) = K,uj\{x). 



(5.10) 
(5.11) 
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From these equations we see that 



uJi = -7^ + coth(x) (5.12) 

is a solution to (|5.11|) and has the correct behavior at infinity. Similar equations for 
the Lie algebra case were found in [T31ll4[ll6H17H26| . This implies the following relation 
for the free energy, 



J_ Fl = -I^ + cothW (5.13) 
1 d dF x, . 

= "2^Mas - coth( 2 ) + coth(x) - 

Here the method of 13] breaks down. In that situation one could trivially integrate to 
get F 1 but here we are unable to. Doing so would amount to writing coth(x) — coth(|) 
as d y^ of some function, something we were unable to do. We will derive Fi using 
the same method as we did in the section 2.2 for the Lie algebra case. One needs to 
find a subleading term in the potential, the term that one has to subtract in order to 
kill the subleading part F\ of the free energy. The first step to this goal is to find the 
leading order resolvent u 



6 Free Energy of Leading Order. 

To find the resolvent ujq(z) in the loop equation it is easier to go back and derive an 
equation of motion. Since we are mostly interested in the Chern-Simons matrix models 
we again restrict our attention to one cut solutions. Let introduce a density function 
Po(u) by, 

[ a ( z — u z + u\ 

Uo(z) = 9s j Po(u) ( coth — h coth — ^— J du. (6.1) 

It is convenient to continue the density function to the negative part of real axis 
p(z) = p(—z). The above definition becomes 



f a z - u 

u (z) = g s J po(u) coth^— du. (6.2) 
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The normalization condition that guarantees the correct behavior of the resolvent at 
infinity is given by 



Po(u)du = 2S. (6.3) 

Lets plug this definition into the loop equation for ujq(z) (|5.1()jl . Subtracting the loop 
equation evaluated at the point z + ie above the cut from the loop equation at the 
point z — ie bellow the cut and taking into account that 



, z-x-k z — x + ie .,. . 
coth coth = 4mo{z — x), (6.4) 



one gets 



2 f a z - u 

—W'(z)=P Po (u)coth——du. (6.5) 

9s J -a " 

The usual way to proceed is to go to a new coordinate U' = e u [3 . In the case of 

SO(2M) Chern-Simons matrix model, the potential is W(z) = z 2 /4 and the equation 
of motion becomes 

--Llog(t/e- 2S ) = P r^jdU', (6.6) 
2g s J e - a U'-U 

where U = e z . Here the normalization condition 



»'} P*Wjfr = 2S (6 - 7) 

has been used. Following |2ll2Zj it is easy to find the function 



<U) = 9s I {^ZJjdU', (6i 



that satisfies the following 



i) vanishes at infinity, 

ii) has a square root cut, 
Hi) v(0 — ie) = 2S, 

iv) v(U - ie) + v(U + ie) = -l/g s log(f/e" 25 ). 
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The only difference from the SU(M) case considered in [3| is that S gets doubled, 

1 + U+ J (1 + U? - AUe 2S 
v(U) = log . (6.9) 

The relationship between this function and the resolvent is 

lu (z) = 2S - 2v(e z ). (6.10) 

Let's discuss the geometry of the Riemann surface given by the resolvent v(u). The 
spectral curve that corresponds to the resolvent is given by 

e v - e~ u + e ~ u ' v+2S - 1 = 0. (6.11) 

Since the resolvent has the property v(u) = v(u + 2iri) the Riemann surface is compact 
in the imaginary direction. The resolvent has the square root cut giving rise to the two 
sheets of the surface. Therefore the Riemann surface looks like two infinite cylinders 
glued together along the cut. The contour around the cut is usually called an A cycle, 
the contour running from a point at infinity on one sheet to a point at infinity on the 
other sheet is called a B cycle. The Riemann surface is depicted in fig. ^ and fig. 121 

From the string theory point of view the curve is part of the mirror B-model ge- 
ometry The branes have disappeared. In terms of type IIA strings there are two 
orientifold 05-planes at points u = izV. These will play a role in the next section 
when we calculate the 0(g s ) part of the free energy. From the string theory point of 
view, the 05-plane contribution to the superpotential is an integral from the location 
of the 05 plane to the point at infinity. The two 05-planes have different D-brane 
charges therefore their contributions are summed with opposite signs. We will see how 
these results emerge from the matrix model. 

Now we are ready to calculate the leading contribution Fq(S) to the free energy. As 
usual it is given by the integral of the resolvent over the B cycle. The integral can be 
expressed in terms of Euler's dilogarithm function (see appendix B for definition and 
useful properties) 

d s F (S) = \ f u (z)dz = - [ v(U)^- = -Li 2 (e~ 2S ) . (6.12) 

^ JB J B U 

Here all infinite and polynomial in S terms are omitted, leaving the worldsheet instan- 
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ton contribution, 



00 -2nS 



71=1 

which agrees with ^H]- The next objective is to calculate F\ and cast it into the form 
of 

7 Free Energy of Order 0(g s ) 

The next objective is to calculate 0(g s ) contribution to the free energy. Although there 
is a derivative relation between ioq(z) and uJ\{z) ()5.12|) . this cannot be integrated to a 
relation between Fq and F\ ()5.13j) . Therefore to find the 0(g s ) part of the free energy 
we use the new method that we have described in the second part of section 2. 



7.1 Fi as a dilogarithm 

To do so we consider the origin of the 0(g s ) term in the free energy. Since the saddle 
point method is used to construct the perturbative expansion in powers of g s , there 
should not be any terms of order of g s unless there is a subleading term in the effective 
matrix model action. To single out such a piece, we add the following subleading term 
to the action, 

Tt5W{u) = 2 !ogsinh 2 u { . (7.1) 

£ 

This is analogous to (j2.11j) in the case of the Lie algebra matrix model. For this new 
potential, W + 5W, we denote the free energy F x . The equation for ujq{z) (|5.1U|) is 
invariant but the equation for uj\{z) ()5.11|) becomes 

}C(z)uji = loq{z)lo\{z) . (7.2) 

Provided that u)i(z) vanishes as z — > ±00 this integral equation has only the trivial 
solution u)i = 0, which leads to Fi = 0. This suggests that ()7.1|) cancels the subleading 
part of the action. So, similar to (J2.19j) . we have, 

- 1 r a 

= -Ft = -Ft + - / p (z) logsinh 2 zdz. (7.3) 
2 Jo 
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In principle, if one knows the density function po(z), the integral can be taken. It is 
more convenient however to have it written as an integral of luq(z), which is known. 
We will see that written in that form, F\ has two different pieces, corresponding to the 
integration along different contours. 
It is easy to see that 

Z Z ~ t - ITT Z7T 

log sinh z = log 2 + log sinh - + log sinh — — . (7.4) 

In what follows we omit all infinite terms and polynomial terms that can be easily 
restored. Similar to [221, we write the logarithms in the integral (j7.3|) as 

z 1 „ f A , x — z , z 



log sinh - = — P / coth dx (7.5) 

2 2 ./q 2 2 



z + in 1 f A x + z z 

logsmh = — / coth ax H — (7.6) 

2 2.4 2 2 



where A is a point at infinity, a UV cutoff. Now one can recognize the resolvent in 
(|7.Hj) . Combining the above and using the fact that po(z) is an even function we get 

Fl ^ S) = ~\{ F l "o( z ) dz + J, uo(z)dz\. (7.7) 

To figure out which branch of the function u)q(z) is to be used in the above integrals or 
in other words on which sheet the point A is located, one has to look at the behavior at 
infinity of the integrals of the coth(^/2) function. The conclusion is that A is a point 
on the physical sheet where ujq(z) vanishes at infinity. The principal value integral is to 
be understood as P J Q = l/2(f Q+ie + f _ ie )- The two contours are drawn on fig. [T] So 
far our discussion in this section has been applicable to one cut solutions of the matrix 
model with arbitrary potentials. From this point we restrict ourselves to quadratic 
potentials only. 

The first integral in ()7.7j) can be written as an integral over the B cycle. To see this, 
we introduce a function y(z) that corresponds to the singular part of the resolvent, so 
the integrals over a cycle of y(z) or of the resolvent are the same. The function y(z) 
has the property of a square root, i.e. it changes the sign when passing over the cut, 
y(z + it) = —y{z — it). The differential y(z)dz is the meromorphic one-form on the 
Riemann surface defined by the resolvent. One has 

u (z) = -2W'(z)+y(z), (7.8) 
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Figure 1 : The physical sheet of the spectral curve. 

-ip £ u Q {z)dz = ~\J "o(z)dz + ~W(A). (7.9) 

To check the result (|7.7jl we calculate Fi(S) for the quadratic potential and compare 
it with the known result. The integral over the B cycle has been taken already in the 
previous section. Therefore 

~P £ u (z)dz = ±U 2 (e~ 2S ) . (7.10) 

*J 

The second integral in ()7.7j) can also be expressed in terms of Euler's dilogarithm 
function (see Appendix B for the details), 

-\£ u (z)dz = ^Li 2 (e" 25 ) + i (Li 2 (e- 5 ) - Li 2 (- e - 5 )) . (7.11) 
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Figure 2: The Riemann surface and and the integration contours that give 0(g s ) part 
of the free energy. 

where again we dropped all polynomial terms. After summing (|7.1U|) and (|7.11|) one 
gets the correct free energy Fi(S) (|4.6j) . 

7.2 Fx as a contour integral 

There is another way of writing (J7.11|) which makes contact with the results of j'i.'jj . 
We want to separate clearly the part coming from 05-planes and the part from the 
expansion of Fo(t). The last one is just the integral over the B cycle. Again it is 
convenient to appeal to the meromorphic one- form y(z)dz 




•A 



1 




A 




) 



uo(z)dz 



2 



y(z)dz + 



y{z)dz 



(7.12) 
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The contour in the second integral is on the second sheet of the Riemann surface, and 
A' is a point at infinity on this sheet (see fig. |2J). Let denote the contour in the first 
integral as O and in the second as 02- From figure El it is clear that O = 0\ + B where 
the contour 0\ goes from in to A'. The conclusion is 

u (z)dz — - ( I (J (z)dz + / u> (z)dz + / u (z)dz j (7-13) 
t 2 \J B J 0l Jq 2 J 

The contour 0\ and the contour O2 run from the positions of two 05-planes to the 

point at infinity which is precisely what is expected from the string theory point of 

view. This consideration so far is valid for arbitrary potentials. If the potential is 

quadratic one has to take integrals of v(u). The change of variables v(u)du = vu'(v)dv 

brings the integrals over 0\ and O2 contours to the same form as in |2*3*j . 

Almost the same calculation can be done for the case of Sp(2M) group. The 

group measure for Sp(2M) has the extra factor Yli sinh 2 Ui which corresponds to the 

additional term ^\ log sinh 2 Ui in the matrix model effective action. To have F\ = 

the term 

Tr5W(w) = -ij^logsinhV (7.14) 

i 

has to be added to the potential. It is clear now that the only difference from SO(2M) 
case is the opposite sign of Fi(S) which agrees with (|4.6|) . 

While it is possible to continue the Chern-Simons partition function to SO(2M +1) 
[TH] and take the large M limit, this partition function does not match the partition 
function of the SO(2M + 1) matrix model. So although the SO(2M + 1) matrix model 
does not corresponds to Chern-Simons theory one can still consider this model. The 
group measure has an extra factor Fj^ sinh 2 y. Repeating the above procedure one gets 

Fi = ^Yl ( lo s 2 ~ y + lo s sinh2 Ut 2 m ~ log sinh2 "2 ) ' ( ^ 7 ' 15 ' ) 

i ^ " ' 

Since the second logarithmic function has changed sign compared to the SO{2M) case 
the two integrals over the B cycle cancel each other and we are left with integrals 
over the contours 0\ and O2. So the part that is proportional to the derivative of Fq 
disappears. This is in contrast with the Lie algebra case, in which the SO(2M + 1) 
and Sp(2M) matrix models have the same free energy. 
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8 N=l SYM In Four Dimensions 



The tree level superpotential is figured out (following [2]) by converting the Haar 
measure into a measure on the Lie Algebra 



J^^smh^-^] ( 2 sink 

i<j ^ ' ^ 



2 / \ 2 

. U{ + M 7 ' 



2 



where, 



\2fe _ 



UK 2 - «?) 2 exp ]T a fc (a k + a+) (it) ) , (8.2) 

i<j \A;=1 



fc(2A;)!' 

i<j v ' 

i?2fc are the Bernoulli numbers. Using the fact that the Newton polynomials 

M 

»=i 

are equal to |TrM fc , where M is an anti-symmetric matrix gauge fixed to the diagonal, 
we see that the four dimensional tree-level superpotential engineered by this construc- 
tion is 



W tree ($) = (8.4) 
8 2 ^ 

k=0 

Only even powers of $ appear, as expected. 

The Chern-Simon's partition function will give the four dimensional low energy 
Wilsonian effective action for Af = 1 SYM with W tTee given by ([8.4)1 . Now in [TH], the 
correct closed string variable was identified as t = g s (2M — 1) (for SO(2M)), which is 
identified with the gluino condensate. So following [TB1IIH! we propose the formula 

BF CS 

W eS = Qds^- + QosG c s - rt, (8.5) 

where Qdq is the total D6-brane charge and Qoe is the total 06-plane charge, r is the 
gauge coupling. Op-plane charge is given by ±2 5_p , -sign for SO, +sign for Sp. We 



^ 2k\ 1 

+ > Tr$ s Tr$ 

^ s AN 



2k-s, 



1 + 
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have introduced Go — aF^ s ' , a is a constant. It is important to use F cs not F MM in 
18.51 because their g s expansions differ, as discussed in the introduction. 

To make contact with results of |2Hj, we look at log(vol(S'0(2M)). This term is 
already within the free energy but we know from [2PJ that it is this term which supplies 
the t log* term to the superpotential. So, expanding in (2M — 1) and keeping only log 
terms, we get 

- log(vol(50(2M)) ~ gfj log t + gj 1 - log t + . . . (8.6) 
so using ()8.5|) . we find 

W eS = ft log* - ±| log* + f ^ - 2FT* ~ rt, (8.7) 
where F pert = F + log(vol(S'0(2M)). So requiring N — 2 vacua, we find that a = 4. 



9 Conclusion 

We have studied the matrix models with Haar measure on SO/Sp in the large M limit, 
for which we have introduced the new form of the resolvent that is compatible with 
the group measure. We have derived a quantum loop equation and for the case of 
quadratic potential, and have found the leading order resolvent. 

We have calculated the 0(g s ) corrections to the 5*0 and Sp Gaussian matrix model 
free energy using a novel method. This method separates in a clear way the leading (of 
order l/g s ) and subleading (of order 0(1)) parts in the effective matrix model action. 
The free energy of the first two orders was expressed as integrals of the leading order 
resolvent over the spectral curve and these integrals were explicitly performed. 

We have found agreement between matrix model and large M Chern-Simons results. 
While the g s expansion of the Chern-Simons theory has a nice worldsheet interpreta- 
tion, which means that the first two orders correspond to sphere and HIP 2 worldsheets, 
the g s expansion of the matrix model free energy mixes the worldsheet contribution 
at each order, essentially due to a shift in the identification of the 't Hooft parameter. 
The derivative relation found in Lie algebra matrix models 

f 1 = (9 1) 

1 4 dS 1 } 
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does not hold but instead we find that 




1 ( dF (S/2 + m) dF (S/2) dF (S) 

2 V dS dS dS 



) 



(9.2) 



Here F 1 contains a contribution from sphere worldsheets as well as HP 2 worldsheets. 

Type IIA string theory on M 3 ' 1 x T*S 3 with the internal geometry orientifolded and 
N D6 branes wrapped on M 3 ' 1 x S 3 engineers an J\f = 1 SO / Sp SYM in four dimensions 
with a certain double trace tree level superpotential which was given. The calculation 
of the leading and subleading free energy in these matrix models or equivalently in 
Chern-Simons theory gives the effective superpotential for this four dimensional SYM. 
This was also discussed. 

Although we have presented the main results with the potential W(z) an arbitrary 
polynomial, the string theory application of these matrix models is known only for 
the quadratic potential. One can convert the matrix model to a Lie algebra matrix 
model with double trace potential and a single trace potential W(Q). This potential 
then corresponds to the tree level superpotential of an J\f = 1 SYM in four dimensions. 
Therefore one knows the four dimensional effective theory but does know the internal 
geometry which constructs this theory. If the potential has higher than quadratic 
powers, the spectral curve will not be a polynomial in e u and e v . Whether or not this 
spectral curve can be related to some B-model geometry is an interesting question to 
address. 

It would be interesting to generalize the impressive work [3U] and solve the SO / 'Sp 
matrix model to all orders by the method of orthogonal polynomials. 
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A Derivation Of the Loop Equation 

In this paper we are just interested in calculating the free energy contribution from 
HP 2 worldsheets. To do this we derive an all worldsheet loop equation. This equation 
will also be derived for Lie group SU (M) for completeness. 
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A.l S0{2M) 

The partition function, after integrating out the off diagonal terms, is given by, 



Z~ f n^nsmh 2 (^^)sinh 2 (^±^)exp (--W(u t )) . (A.f) 
•* i=i j^i ~ ~ \ 9s J 

We perform the infmitessimal change of co-ordintes, 

Ui^Ui + e (^coth(^p) - coth(^p)^ . (A.2) 

and demand that the partition function is invariant under this transformation. At first 
order this yields the following constraint, 



(^cosech 2 (^-y^) + cosech 2 (^4^)^ 
+ E (-th(^p) - coth(^)) (coth(^) + coth(*±^)) 
2 / coth(^) - coth(^)) W'(u0\ = 0. (A.3) 



9s 

Now we need the 2 identities, 



^ cosh(^) 



^sinh(^)sinh(^) 

M 

cosech( )cosech( -) + 2 > coth( ) coth( -) 

Z Z Z Z 

1=1 i^j 

M 

= Y: coth(^p) coth(^) - M 2 (A.4) 



and 
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M 

E 



cosh(^ 



; sinh(^) sinh(^) 



1,3- 



M 



cosech( ^ 

i=l 
M 



xosech( — — — 



2 coth( 



)coth( J ) 



— ;COth(^)-M 2 . 



With a little more work, one can also show that 



(A.5) 



-u 2 (x) — # s coth(x)o;(x) + g s S 



2 ^ 



M r 



+£ 2 £J ( COth( 



( , )t]l (^±^)coth(^^ 



,X — U 



2 



coth( ) coth( 



.T — Mj 



coth( 



I) ) C oth(^r^) + S 2 - ^M. (A.6) 



Now we multiply ()A.3|) through by g 2 and employ (|A.41 IA.51 IA.6|) to get the final 
form of the loop equation, 



\uj\x) - g s coth(;rV0r) + 2g s S - 2S 2 - Ku(x) + ^-^-^(x) = 0, (A.7) 



where the linear operator /C acts as 



e ti C °^ )W ' {z)f{z) - 



(A.i 



A. 2 SU(M) 

Just as above, we integrate out the off diagonal components and demand reparametri- 
sation invariance under 

Ui — > Ui + ecoth(— -— (A. 9) 
Since the derivation is much simpler than for SO(2M), we just state the result, 



^(x)-2^-^) + f^(x) = 0. 



(A.10) 
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B Dilogarithm Identities 



Euler's dilogarithm function Li 2 is defined as the integral 



U 2 (z) = f'^f^dt, (B.1) 



or as the power series 



00 u 

z K 



71=1 

Among the many amazing properties of this function, we will use the following (for 
a review see for example [3*T] ) 



Li 2 (^) + Li 2 (-^) = iLi 2 (^ 2 ), (B.3) 

Li 2 (l -z) + Li 2 (l - z- 1 ) = l - log 2 z, (B.4) 

U 2 (z) + Li 2 (l - z ) = ^r- (log*)(log(l - z)). (B.5) 

o 



The integral from m to A can be taken, with the result 



If 1/1 

- J u (z)dz = -- log 2 (~e- s ) - log (-e 5 ) log (l + e" 5 ) 

+Li 2 (e~ 2S ) - Li 2 (-e- s ) - Li 2 (l + e s ) ) . (B.6) 



To cancel the product of logarithms one has to apply (|B.4J) to Li 2 (l + e s ) to get 
Li 2 (l + e _s ), then using (|B.5|) convert it to Li 2 (— e _s ). After this manipulation and 
omitting all polynomial terms in S, one has 

i J A u (z)dz = -1 (-2Li 2 (-e- s ) + Li 2 ( e - 25 )) . (B.7) 
The last step is to apply the property (jB.3|) and get (|7.11|) . 
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